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The  satellite  1971-54A  entered  a  near-circular  orbit  with  period  95.9  minutes 
and  inclination  90.2°.  Between  1972  and  1978  the  orbit  passed  slowly  through 
15th-order  resonance,  when  the  track  over  the  Earth  repeats  after  15  revolutions, 
and  the  15th-  and  30th~order  harmonics  in  the  geopotential  may  produce  substantial 
orbital  perturbations. 


The  values  of  orbital  inclination  and  eccentricity  from  269  weekly  US  Navy 
orbits  between  November  1972  and  January  1978  have  been  analysed  to  determine  12 
lumped  harmonic  coefficients  of  order  15  and  30.  The  analysis  of  inclination 
yields  15th-order  coefficients  accurate  to  1.5%  and  2.8%,  and  30th-order  coef¬ 
ficients  accurate  to  7%.  The  analysis  of  eccentricity  gives  two  15th-order 
coefficients  accurate  to  3%  and  4%. 


These  lumped  harmonic  coefficients  are  used  to  test  the  accuracy  of  the 
Goddard  Earth  Model  10B,  which  is  complete  to  order  and  degree  36.  The  agreement 
with  GEM  10B  is  excellent,  for  both  15th  and  30th  order,  and  shows  that  GEM  10B 
is  more  accurate  than  was  expected. 

The  12  values  of  lumped  harmonics  obtained  give  12  linear  equations  between 
individual  coefficients  of  order  15  and  30,  which  will  be  used  in  a  future  solu¬ 
tion  for  the  individual  coefficients,  c 
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INTRODUCTION 


T!ie  Earth's  gravitational  potential  is  usually  expressed  as  a  double 
infinite  series  of  tesseral  harmonics  of  degree  l  and  order  m  (41)  .  The 
I’rdc”  m  of  the  harmonic  specifies  variations  between  one  meridian  and  another, 
and  if  the  Earth  were  sliced  through  the  equatorial  plane,  a  harmonic  of  order  m 
would  express  itself  as  a  shape  with  m  oscillations  as  the  longitude  goes 
through  360°:  thus  there  would  be  'humps'  at  intervals  of  360°/m. 

The  orbit  of  an  Earth  satellite  exhibits  nth-order  resonance  when  its 
or  >L»al  period  is  such  that  its  path  over  the  Earth  repeats  after  m  revolutions. 
Successive  tracks  over  the  Earth  are  then  spaced  at  360°/m  in  longitude,  and  the 
harmonics  of  order  m,  2m,  3m,  ...  have  the  same  effect  on  the  orbit  on  succes¬ 
sive  transits.  So  their  effects  build  up. 

In  practice,  a  satellite  close  to  the  Earth  suffers  appreciable  air  drag, 
as  a  result  of  which  its  orbital  period  steadily  decreases,  and  the  satellite  is 
slowly  drawn  through  the  resonance  by  drag.  The  effects  of  the  resonance  may  last 
for  several  years,  however,  for  a  low-drag  satellite,  and  the  orbit  may  be 
altered  by  as  much  as  5  km. 

A  thorough  analysis  of  such  changes  can  give  accurate  values  for  'lumped 
harmonics'  of  order  m  ,  a  lumped  harmonic  being  the  linear  sum  of  a  series  of 
individual  harmonic  coefficients  of  order  m  and  degree  5.^,  1^+2,  +  4,  ..., 

where  =  m  or  m  +  1  .  The  analysis  also  usually  yields  values  of  lumped 
harmonics  of  order  2m,  of  poorer  but  still  useful  accuracy;  and  sometimes  gives 
lumped  harmonics  of  order  3m.  In  these  lumped  harmonics  there  are  often  signi¬ 
ficant  contributions  from  individual  harmonics  up  to  degree  50  or  more. 

Accurate  values  of  lumped  harmonics  are  now  particularly  useful,  because 

they  provide  an  exacting  independent  test  of  the  accuracy  of  the  comprehensive 

models  of  the  gravity  field  determined  from  other  data.  Until  recently,  the 

tests  were  often  inconclusive  because  the  models  did  not  go  beyond  order  and 

degree  36  -  the  limit  of  GEM  10b',  which  was  derived  from  some  800000  accurate 

observations  of  Earth  satellites,  including  about  200000  each  of  laser,  Doppler 

and  camera  observations,  together  with  5°  x  5°  surface  gravity  measurements  and 

radio  altimeter  data  on  700  passes  of  the  Geos  3  satellite.  By  the  addition  of 

28000  1°  x  1°  areas  of  sea-surface  geoid  derived  from  2300  passes  of  altimetry 

2 

from  Geos  3,  a  new  geopotential  model,  GEM  IOC,  has  been  derived  ,  which  goes  up 
to  order  and  degree  180,  but  is  the  same  as  GEM  10B  up  to  order  and  degree  36. 

The  existence  of  GEM  IOC  is  a  spur  to  new  and  more  accurate  determinations  of 
lumped  harmonics. 
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Many  examples  of  orbital  resonance  have  been  analysed  in  recent  years. 

3  4  5 

Individual  harmonics  of  15th  order  ’  and  14th  order  ,  and  also  tentative  values 
of  30th-order  coefficients  ,  have  been  determined  from  analysis  of  a  number  of 
orbits  at  different  inclinations.  In  addition,  lumped  harmonics  have  been 

determined  from  examples  of  29:2  resonance^  (29  revolutions  in  2  days),  and 

89io  1112 

31:2  resonance  1  ’  ,  as  well  as  from  11th-  and  13th-order  resonances  ’  .  Many 

studies  of  near-resonant  orbits  have  also  been  made. 

The  most  accurate  values  of  lumped  harmonics  are  likely  to  be  obtained  from 
satellites  of  low  drag,  and  1971-54A  was  chosen  for  analysis  because  it  was 
appreciably  affected  by  !5th-order  resonance  for  more  than  5  years,  from  mid- 1972 
to  the  end  of  1977. 

Through  the  courtesy  of  the  US  Naval  Resea^'h  Laboratory,  we  have  received 
US  Navy  orbital  elements  for  1971-54A,  and  the  Navspasur  observations  of  the 
satellite,  at  weekly  intervals  since  launch.  The  present  analysis  is  entirely 
based  on  269  weekly  US  Navy  orbits,  for  epochs  between  12  November  1972  and 
14  January  1978.  The  orbital  elements  used  are  the  inclination  and  eccentricity. 

As  it  happens,  high-degree  harmonics  have  relatively  little  effect  on 
I971-54A,  because  the  orbit  is  polar,  so  the  results  provide  a  test  for  GEM  1  OB 
rather  than  GEM  10C. 

2  THE  SATELLITE  1 97 1 -54A  (SESP-1)  AND  ITS  ORBIT 

The  satellite  designated  1971-54A  was  launched  by  the  USAF  on  8  June  1971 

into  a  near-circular  orbit  of  inclination  90.2°,  with  perigee  and  apogee  heights 

near  550  and  590  km.  The  satellite  consisted  of  a  Burner  2  rocket  with  a  payload 

attached:  called  the  'Space  Experiments  Support  Program  1'  (SESP-1),  its  purpose 

1  3 

is  said  to  have  been  'to  test  an  infra-red  celestial  mapping  sensor  system' 

The  size,  shape  and  mass  of  the  satellite  are  not  exactly  known,  but  it  was 
1 4 

probably  about  3  m  long  and  1.3  m  in  diameter,  with  a  mass  of  about  260  kg. 

The  satellite  was  launched  at  a  time  when  the  solar  activity,  and  hence  the  upper- 
atmosphere  density,  was  decreasing  after  the  maximum  of  1968-70.  Consequently, 
the  rate  of  orbital  decay  was  very  slow  between  1972  and  1977,  while  the  solar 
activity  remained  low.  The  decrease  in  orbital  period  between  1971  and  1979  is 
shown  in  Fig  I  . 

By  great  good  fortune,  the  15th-order  resonance  occurred  during  the  time  of 
lowest  drag,  thus  maximizing  its  effects.  The  orbit  has  been  analysed  between 
12  November  1972  and  14  January  1978,  as  indicated  in  Fig  1:  throughout  this  time 
the  effects  of  !5th-order  resonance  were  appreciable  because  the  anomalistic 
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period  was  within  0.1  minute  of  its  resonant  value,  95.68  minutes.  The  analysis 
was  shopped  at  14  January  1978  because  of  the  great  increase  in  drag  as  a  result 
of  the  increased  solar  activity.  It  will  be  seen  from  Fig  1  that  for  a  short 
time  early  in  1975  there  was  an  increase  in  orbital  period  because  the  effect 
of  the  resonance  was,  temporarily,  strong  enough  to  overcome  the  effect  of  drag. 
Although  I971-54A  was  too  low  in  the  atmosphere  to  become  trapped  in  resonance, 
it  hesitated  before  escaping,  and  nearly  turned  back  to  experience  resonance  for 
a  second  time. 

3  PREPARATION  OF  THE  ORBITS  FOR  THE  RESONANCE  ANALYSIS 

The  resonance  analysis  is  conducted  with  the  aid  of  the  THROE  computer 
program*  ■*,  which  requires  orbital  elements  in  the  5-card  format  of  the  PROP 
orbit  refinement  program1^.  The  US  Navy  orbits  differ  from  PROP  orbits  in  two 
ways : 

(a)  The  epoch  is  at  the  last  ascending  node  before  midnight,  rather 
than  at  midnight  as  in  PROP; 

(b)  The  values  of  eccentricity  have  a  part  of  the  odd  zonal  harmonic 
oscillation  removed,  and  the  argument  of  perigee  also  consequently 
suffers  modification. 

The  269  US  Navy  orbits  were  converted  to  PROP  elements  at  the  nearest  midnight 
using  the  program  ELTRAN2,  which  also  restores  the  perturbation  in  eccentricity 

on  the  assumption  that  the  values  of  the  geopotential  coefficients  J  and  J,. 

"”6  “6  ^ 
used  were  =  -2.551  *  10  and  Jr  =  -0.191  *10  •  After  conversion  by 

ELTRAN2,  the  269  values  of  eccentricity  e  and  argument  of  perigee  co  closely 

followed  the  expected  basic  circular  pattern  of  variation  -  that  is,  a  circle 

in  the  (e  cos  u>  ,  e  sin  w)  cartesian  plane.  So  the  conversion  was  considered 

successful . 

Before  the  resonance  can  be  analysed,  other  perturbations  must  be  removed 

from  the  orbital  elements.  For  the  inclination,  the  most  important  perturbations 

are  those  due  to  lunisolar  gravitational  attraction  and  the  zonal  harmonics  in 

the  geopotential.  These  were  calculated  by  running  the  PROD  computer  program'^ 

with  I -day  integration  steps  through  the  1889  days  of  the  analysis,  with  restarts 

at  intervals  of  about  100  days,  because  the  running  integration  gradually  drifts 

off  as  a  result  of  the  air-drag  and  resonance  effects  not  included  in  PROD.  The 

1 8 

lunisolar  gravitational  perturbations  are  themselves  resonant  for  a  polar  orbit  , 
and  therefore  near-resonant  for  1971-54A  at  an  inclination  of  90.2°.  During  the 
5  years  of  the  integration  by  PROD  the  inclination  irregularly  but  inexorably 
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increased  by  about  0.05°.  Since  this  perturbation  was  so  large,  its  zero  was 
taken  to  occur  at  a  time  close  to  resonance,  so  that  the  maximum  departure  of  the 
adjusted  values  from  the  observed  values  would  be  as  small  as  possible.  To 
remove  the  effects  of  lunisolar  and  zonal  harmonic  perturbations,  the  amount  to 
be  added  to  the  observed  value  of  inclination  was  0.0170°  initially;  0.0091°  at 
MJD  42000;  -0.0078°  at  MJD  42500;  -0.0143°  at  MJD  43000;  and  -0.0296°  at  the  end 
of  the  analysis. 

The  second  important  perturbation  of  the  inclination  is  that  due  to  atmos¬ 
pheric  rotation.  For  a  circular  polar  orbit  with  a  period  decreasing  by  AT 
minutes,  the  expected  decrease  in  inclination  due  to  an  atmosphere  rotating  at 

A  rev /day  is  0.007A  AT  degrees.  Here  AT  =  0.2  over  the  5  years,  so  the  inclina- 

.  .  1 9 

tion  decreases  by  0.00I4A  degrees;  previous  results  indicate  an  average  atmos¬ 
pheric  rotation  rate  of  0.8  +  0.1  rev/day  at  heights  near  450  km,  so  the  expected 
decrease  in  inclination  is  0.0011°  +  0.0001°  over  the  5  years.  This  perturbation 
is  computed  within  THROE,  together  with  that  due  to  the  precession  of  the  Earth's 
axis.  Perturbations  due  to  Earth  and  ocean  tides  are  expected  to  be  about  50  m 
at  most,  and  need  not  be  considered,  since  the  ar  -uracy  of  the  observational 

values  of  inclination  is  about  200  m.  Perturbations  due  to  solar  radiation 

.  2 
pressure  for  1971-54A,  which  has  a  mass/area  ratio  of  about  200  kg/m  ,  are 

estimated  at  0.00004°,  and  are  also  negligible. 

By  far  the  most  important  perturbation  in  eccentricity  is  that  due  to  the 
odd  zonal  harmonics  in  the  Earth's  gravitational  potential,  which  produce  an 
oscillation  with  the  same  period  as  w  (about  100  days)  and  an  amplitude  of 
0.0014,  which  is  equivalent  to  about  10  km  in  perigee  distance.  This  perturba¬ 
tion  is  removed  within  THROE  using  a  20-coefficient  set  of  zonal  harmonics,  the 
odd  zonal  harmonics  being  those  of  Ref  20.  Since  the  eccentricity  is  small 
(its  average  value  being  0.0022)  the  lunisolar  perturbations  to  e  should  also 
be  small.  They  were  evaluated  by  running  PROD  with  and  without  lunisolar  per¬ 
turbations:  it  was  found  that  the  maximum  difference  in  e  between  the  two  runs 
was  0.000001,  which  is  negligible.  Earth  tide  and  radiation  pressure  perturba¬ 
tions  to  e  should  also  be  negligible. 

If  the  perturbations  in  inclination  and  eccentricity  are  successfully 
removed,  the  remaining  variations  in  these  two  orbital  parameters  should  be 
caused  only  by  the  effects  of  15th-order  resonance. 
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THEORY  FOR  GENERAL  6:a  RESONANCE 


The  longitude-dependent  part  of  the  Earth's  gravitational  potential  at  an 

2 1 

exterior  point  (r,0,A)  may  be  written  in  normalized  form  as 


V 
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i  ( 

ZZ($ 


j.x  i  _  m  f 
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1=2  m= I 


0)1 


C.  cos  mA 
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(1) 


where  r  is  the  distance  from  the  Earth's  centre,  9  is  co-latitude,  X  is 

longitude  (positive  to  the  east),  u  is  the  gravitational  constant  for  the  Earth 

(398600  km  /s"),  R  is  the  Earth's  equatorial  radius  (6378.1  km),  P™(cos  8)  is 

the  associated  Legendre  function  of  order  m  and  degree  Jl  ,  and  C.  and  S„ 

Hm  Jim 

are  the  normalized  tesseral  harmonic  coefficients.  The  normalizing  factor  N„ 

2  1 

is  given  by 

„2  _  2(21 1  +  DU  "  m)! 

Njim  -  - rmsyi —  •  (2) 


When  a  satellite  orbit  experiences  B:a  resonance  -  that  is,  when  the 
satellite  makes  6  revolutions  while  the  Earth  makes  a  revolutions  relative  to 
the  satellite's  orbital  plane  -  the  rates  of  change  of  inclination  and  eccen¬ 
tricity  depend  on  the  resonance  angle  $  ,  defined  by 

t  =  a(w  +  M)  +  8(J1  -  v)  ,  (3) 


where  M  is  the  mean  anomaly,  SI  is  the  right  ascension  of  the  ascending  node, 

and  v  is  the  sidereal  angle.  The  perturbation  produced  in  i  and  e  by  a 

.  .  -  -  22  23 

relevant  pair  of  harmonic  coefficients,  C.  and  S„  ,  may  be  written  ’ 
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Here  F_  is  Allan's  normalized  inclination  function  ,  G„  is  a  function 
Imp  Jlpq 

of  eccentricity  e  for  which  explicit  forms  have  been  derived  by  Gooding  ,  ^ 
denotes  'real  part  of'  and  j  =  /-T  .  The  indices  y,  p,  k  and  q  are  integers, 
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with  y  taking  the  values  1,  2,  3,...,  and  q  the  values  0,  ±1,  ±2,...  .  The 
equations  linking  l,  m,  k  and  p  are:  m  =  yB;  k  =  yet  -  q;  2p  =  £  -  k  . 

For  a  particular  0:a  resonance,  the  relevant  value  of  m  in  the  C.  , 

_  x-m 

S„  coefficients  is  decided  by  the  choice  of  y  ,  since  m  =  yS>  .  There  is  a 

x.m 

series  of  possible  values  of  l  -  any  values  such  that  £  i  m  and  (£  -  k)  is 
even.  The  successive  C  ,  coefficients  that  arise  (for  given  y  and  q) 

may  usefully  be  combined  in  'lumped'  form  and  written  as 


q?’kc, 

^£  £m 


-  !>' 


where  £  increases  in  steps  of  2  from  its  minimum  permissible  value  £  ;  the 

q  k  0 

are  functions  of  inclination  that  can  be  taken  as  constant  for  a  particular 
satellite;  and  Q^’k  =  1  when  £  =  £^  . 

5  ANALYSIS  OF  THE  RESONANT  VARIATION  OF  INCLINATION  FOR  I971-54A 

5. 1  The  form  of  the  equation  for  di/dt 

.  ( 1  £e)  I  **  I 

In  equation  (4)  for  di/dt  ,  the  functions  G.  are  of  order  A2? — 

*pq  ( I  q I ) ! 

and,  since  e  is  so  small  (^0.002)  for  I971-54A,  the  q  =  0  terms  in  equation  (4) 
should  be  the  largest.  For  15th-order  resonance,  with  8  =  15  and  a  =  1,  and 
with  only  the  (y,q)  =  (1,0), (2,0)  and  (3,0)  terms  written  explicitly,  equation  (4) 
becomes : 


di  _  n  / R^j 
dt  sin  i  \aj 


(15  -  cos  i 


l)|^!5,l5,7jC15  Sln  $ 

*  2(^5f30.30,l4|5° 

*  3(!)30f45,«,2l{E4 


C^’1  sin  $  -  s°’!  cos  $ 


sin  2$  -  S°3’Q2  cos  2$ 


C°l3  sin  3$  -  S°’3  cos  3$ 
45  45 


( ^ He) ' ^  cos  ,  ^  , 

+  terms  in  — Vi — i-tt"  •  (y$  ~  qw) 

(  q  ) !  sm 


The  coefficients  F^  ^  F^q  ^  etc,  and  the  multiplying  factors 
Q^’  ,  have  been  calculated  with  the  aid  of  the  computer  program  PROF  for  1971-54A, 

taking  mean  values  of  the  elements  at  a  time  near  the  exact  resonance,  namely 
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a  =  6930  km,  e  =  0.0022  and  i  =  90.21°.  The  numerical  values  of  the  can 

_q,k  l 

be  used  to  make  a  rough  estimate  of  the  order  of  magnitude  of  C  in  the  vari- 

m 

ous  (Y.q)  terms  in  equation  (7),  on  the  assumption  that  the  individual  coef- 

-  -52  -9>k 

frcients  C.  ,  S„  are  of  order  10  /I  .  We  assume  C  is  of  order 
Urn  £m  m 

10-%  (Q’’k/t2)2'! 


By  inserting  numerical  values  for  the  F  functions  and 


R/a  (=  0.9204),  the  relative  magnitudes  of  thevarious  (y,q)  terms  in  equation  (7) 
can  be  calculated.  The  results  are  shown  in  Table  1,  the  leading  term  being 
scaled  up  to  100:  the  explicit  forms  of  the  (y,q)  =  (1,±1)  terms  were  taken  from 
Ref  24. 


Table  1 


Relative  orders  of  magnitude  of  (Y,q)  terms  in  (7),  for  I971-54A 


l 

2 

3 

0 

100 

10 

2 

i 

0.  1 

- 

- 

-l 

1 

— 

“ 

Terms  marked  with  a  dash  and  those  with  jq|  2  or  Y  5-  4  are  not  considered 
large  enough  to  be  worth  evaluating. 

Table  1  shows  that  the  three  terms  given  explicitly  in  equation  (7)  should 
be  the  largest.  If  the  main  term  can  be  evaluated  with  an  accuracy  of  2%  or 
better,  it  may  be  necessary  to  include  the  (Y>q)  =  (3,0)  term  in  the  fitting. 

5 . 2  The  resonance  angle  $ 

From  equation  (3),  the  resonance  angle  <i>  for  15th-order  resonance  is 
given  by 

$  =  u  +  M  +  I5(«  -  v)  .  (8) 

The  variations  of  $  and  4>  during  the  5  years  of  the  analysis  are  shown  in 
Fig  2. 

The  variation  of  $  follows  the  usual  near-parabolic  form,  but  the 
'hesitation'  after  resonance  is  clearly  visible  in  Fig  2. 


I 


For  a  satellite  of  fairly  high  drag,  the  increase  of  $  is  usually 

nearly  linear,  with  variations  in  slope  in  response  to  variations  in  drag.  Here, 

however,  4>  pursues  a  sinuous  course  through  the  resonance,  with  its  slope 

dependent  on  <J>  as  well  as  on  the  drag.  The  equation  for  the  variation  of  M 
22 

due  to  resonance  shows  that  M  -  and  hence  <J>  -  tends  to  decrease  whenever 
di/dt  is  negative.  The  numerical  values  of  C^'  and  *  for  I971-54A  are 

such  that  di/dt  is  negative  wherever  -18°  <  $  <  162°  ,  and  has  its  maximum 

negative  value  when  $  =  72°.  Inspection  of  the  values  of  $  in  Fig  2  shows 

that  the  resonance  effect  is  reducing  M  ,  and  hence  &  ,  for  530  days  between 

MJD  42054  and  42583,  when  57°  <  $  <  162°  .  It  was  only  at  the  end  of  this 
time  interval,  in  June-July  !975,  when  air  density  fell  to  its  minimum  value  in 
the  solar  cycle,  that  the  effects  of  resonance  were  able  to  overcome  drag.  If 
the  mininum  density  had  occurred  a  little  earlier,  the  resonance  effect  could 
have  pulled  the  satellite  back  to  resonance,  where  it  might  have  become  trapped 
until  the  density  increased  as  the  solar  activity  began  its  increase  in 
September  1977. 

The  points  where  the  maximum  (negative)  effect  of  resonance  on  4>  is 
exDected  -  ie  $  =  432°,  792°, 1 1 52°,  .  . .  -  are  marked  on  Fig  2,  and  it  can  be 
seen  that  they  correspond  to  the  minima  in  the  slope  of  $  . 

If  the  effect  of  drag  could  be  independently  assessed  and  removed,  the 
variation  of  i  would  provide  an  easy  and  accurate  method  of  determining  lumped 
harmonics:  but  there  seems  little  prospect  of  being  able  to  estimate  the  drag 
with  adequate  accuracy  by  other  methods. 

The  other  obvious  conclusion  to  be  drawn  from  Fig  2  is  that  the  evaluation 
of  air  density  from  a  satellite's  orbital  decay  rate  should  not  be  attempted 
near  a  resonance,  unless  the  effect  of  resonance  is  very  much  smaller  than  that 
of  drag. 

5 . 3  The  fitting  of  the  variations 

The  analysis  of  the  variation  in  inclination  presented  no  problems,  and 
Fig  3  shows  the  theoretical  curve  fitted  by  THROE  to  the  269  values  after 
removal  of  the  four  perturbations  specified  in  section  3  (lunisolar,  ronal 
harmonic,  atmospheric  rotation  and  precession  of  the  Earth's  axis).  The  curve 
fitted  is  that  specified  by  equation  (7)  with  (y,q)  =  (1,0), (2,0)  and  (3,0). 

The  a  priori  error  assigned  to  the  values  of  inclination  was  0.003°,  but  the 

residuals  are  nearly  all  less  than  0.003°,  so  that  the  weighted  residuals  are 

2 

less  than  I  and  e  =  0.55  (where  e  is  the  sum  of  the  squares  of  the  weighted 
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residuals,  divided  by  Che  number  of  degrees  of  freedom).  Thus  the  a  posteriori 
estimate  of  accuracy  is  0.0016°,  equivalent  to  200  m  in  position. 

The  fitting  of  the  theoretical  curve  in  Fig  3  is  as  nearly  perfect  as 
is  ever  seen  in  the  real  world.  It  is  particularly  significant  that  the  ampli¬ 
tudes  of  the  minor  oscillations  before  and  after  resonance  are  correctly  modelled 
by  the  theoretical  curve,  thus  banishing  any  doubts  about  the  theory  which  may 

have  arisen  from  poor  modelling  of  the  subsidiary  oscillations  in  some  previous 
25 

analyses  .  The  phase  of  the  minor  oscillations  after  resonance  lags  behind 
Lli.it  of  the  theoretical  curve,  but  this  can  be  partially  explained  by  the  fact 
that  the  epoch  of  the  US  Navy  orbits  is  al  Lhe  end  of  the  7  days  of  observations, 
so  that  the  inclination  values  might  arguably  be  regarded  as  applying  to  an 
epoch  days  earlier. 

The  main  decrease  in  inclination  in  Fig  3,  from  90.25°  in  November  1973  to 
90.15°  in  October  1975,  is  equivalent  to  a  positional  change  of  12  km. 

5 . 4  Values  of  the  lumped  coefficients 

-0,1  _0,1 

The  fitting  shown  in  Fig  3  gives  numerical  values  of  and 

having  standard  deviations  of  27.  and  3%  respectively;  so,  on  the  basis  of  Table  1  , 
it  is  questionable  whether  the  45th~order  terms  are  needed.  A  fitting  with 
(y,q)  =  (1,0)  and  (2,0)  only  was  therefore  made:  the  value  of  e  was  0.552, 
almost  the  same;  nearly  all  the  weighted  residuals  were  within  +0.1;  and  the 
standard  deviations  of  the  15th-  and  30th-order  coefficients  were  smaller. 

Table  2  gives  the  values  of  the  lumped  coefficients  obtained  from  the  two 
fittings. 

Table  2 


Values  of  lumped  harmonics  obtained  from  fitting  of  i 


1 

- — - — T 

(r.q) 

pairs  included 

| 

9  0,1 

10  5,5 

q  o,: 
l0C)0 

9  0.2 
,0S')0 

9  0,3 
,0C45 

9,0,3 

10  S45 

0,0) 

(2,0)  (3,0) 

-16.3  :  0.3) 

-3.2  t  0. 18 

-f, .  7  •  : . ; 

1 !  .  7  ?  !  .0 

6.0  ♦  3.4 

-2.5  i  3.0 

0,0) 

(2,0) 

—  1 6 .  ■»  t  0 .  J  •  * 

-5.4  !  0.IS 

-8.2  t  0.6 

1  1 . 1  ♦  0.8 

- 

- 

(The  standard  deviations  of  the  15th-order  coefficients  are  given  to  two  decimal 
places  so  as  to  allow  comparisons  between  them.) 


Since  the  45th-order  terms  are  ill-determined  and  the  (1,0), (2,0)  solution 

gives  better  standard  deviations,  it  is  probably  to  be  preferred.  However,  the 
-0  3 

C  ’  term,  when  multiplied  by  the  ratio  of  the  F  terms  etc  in  (7)  to  provide 

.  .-01 

a  measure  of  its  effect  relative  to  C  ’  ,  has  a  value  1.1  ±  0.6,  which  exceeds 

0  1  ^ 

the  standard  deviation  of  Cj’  :  so  it  could  be  argued  that  the  45th-order  terms 
are  needed  in  the  fitting.  Thus  it  seemed  best  to  give  both  solutions:  it  is  very 
satisfactory  that  they  agree  to  within  the  sum  of  their  standard  deviations. 

The  values  of  both  the  15th-  and  30th-order  coefficients  in  the  second  solution 
are  nominally  the  most  accurate  yet  obtained  from  any  satellite. 

Now  that  the  actual  values  of  the  lumped  harmonics  are  known,  the  predic¬ 
tion  of  their  orders  of  magnitude  in  Table  1  can  be  checked.  The  actual  ratios 
of  the  ( 1 ,0) : (2,0) : (3,0)  terms  are  31:11:2  rather  than  100:10:2  in  the  order- 

of  magnitude  prediction  in  Table  I.  So  the  15th-order  terms  are  relatively 

-5  2 

smaller  than  expected  on  the  simple  10  /£  prediction. 

6  ANALYSIS  OF  THE  RESONANT  VARIATION  OF  ECCENTRICITY  FOR  1971-54A 

6 .  1  The  form  of  the  equation  for  de/dt 

In  analysing  the  eccentricity,  we  have  first  to  decide  which  of  the  (y,q) 
terms  in  equation  (5)  are  likely  to  be  important.  Assuming  that  G  ^  is 

order  (|  le)  ^  q  ^/(  |  q  [ ) ! ,  the  term  G£pq|~ - — 2^— j  in  (5)  is  of  order  £ye 


when  q  =  0;  of  order  when  |q|  =  1  ;  of  order  e  when  |q|  =  2  ;  and  of 

order  j(Ji)^e^  '^^(|q  -  1  | )  !  for  general  |  q  |  >  1  .  Since  the  mean  value 
of  e  for  1971-54A  is  0.0022,  the  |q|  =  1  terras  are  expected  to  be  dominant. 
With  the  same  assumptions  as  in  section  5.1,  the  relative  magnitudes  of  the 
relevant  (y,q)  terms  for  1971  — 5 4A  are  given  in  Table  3,  the  explicit  algebraic 
forms  for  the  terms  again  being  taken  from  Ref  24. 


Terms  marked  with  a  dash  and  those  with  |q|  3  or  y  4  ,  are  not  considered 

large  enough  to  be  worth  evaluating. 


The  very  small  value  of  the  (y,q)  =  (1,1)  term  in  Table  3  is  unusual, 
and  arises  because  its  multiplying  factor  F  ,  0  has  a  zero  at  i  =  90°  and 

q  I  O  j 

is  therefore  very  small  at  i  =  90.2  .  The  disparity  between  the  (1,1)  and  (1,-1) 
terms  does  lead  to  some  problems  in  the  fitting. 


From  Table  3  it  appears  that  the  only  terms  needed  in  the  fitting  will 
probably  be  those  with  (y,q)  =  (1,±1)  and  (2,+l),  and,  if  so,  equation  (5)  for 
de/dt  becomes 


de 

dt 


n(f)  ['  Tf16,15,8{§'l5°  Sin($  '  U)  +  5!s°  COS(*  "  4 

+  -£*16,15, 7K5’2sin(*  +  “)+  5i5,2c0s($  +  4 
-  17(!j5*3.,30,.5{5io'  Sin(2$  -  U)  +  C3o’  C°S(2$  ~  4 

+  *31 ,30, 1 4 ^^30 * 3  Sin(2*  +  W)  +  C~30’3  C°S(2$  +  4]  •  (9) 


This  will  be  the  basic  equation  for  fitting  the  observed  values,  though  other 
terms  will  be  added  and  dropped  on  a  trial  basis. 

6.2  The  fitting  of  the  variations 

The  269  values  of  eccentricity  were  fitted  with  equation  (9)  using  THROE, 
zonal  harmonic  and  atmospheric  perturbations  being  removed  within  THROE.  The 
fitting  was  unacceptable:  the  maximum  weighted  residual  was  8.4;  there  were 
several  long  runs  of  residuals  of  the  same  sign;  and  the  overall  measure  of  fit, 
c  ,  was  3.3.  Re-runs  with  (y,q)  =  (3,±1)  terms  added,  and  then  with  ( 1 , ±2 )  terms 
added,  somewhat  reduced  e  ,  to  between  2.5  and  2.8,  but  the  fitting  was  still 
unsatisfactory.  Finally,  (y,q)  =  (0,1)  terms  were  added  to  the  basic  four  pairs: 
this  introduces  terms  (A  sin  u  -  B  cos  w) ,  which  do  not  arise  in  the  resonance 
equation  (5),  and  their  inclusion  represents  an  attempt  to  remove  any  spurious 
oscillations  of  the  same  frequency  as  u>  and  unconnected  with  the  resonance. 

The  (y,q)  =  (0,1)  terms  added  were  well  determined,  e  was  reduced  from  3.3  to 
2.2,  and  the  fitting  was  much  more  satisfactory.  So  the  run  with  (y,q)  =  (0,1) 
terms  added  can  be  regarded  as  "the  best  of  a  bad  bunch",  and  on  that  basis  it 
is  shown  in  Fig  4. 
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The  fitting  of  the  points  in  Fig  4  is  good  in  one  respect:  the  detail  in  the 
observational  values  is  well  mirrored  by  the  curve,  which  has  ups  and  downs  in 
nearly  perfect  unison  with  the  data.  But  the  curve  is  usually  either  above  or 
below  the  general  run  of  the  values  by  a  margin  too  large  to  be  acceptable. 

If  we  ignore  its  defects  for  the  moment.  Fig  4  has  some  interesting 
features.  First,  it  shows  that  the  variation  in  eccentricity  near  the  <t  =  to 
resonance  is  much  smaller  than  the  variation  near  the  d>  =  -  to  resonance:  this 
is  as  expected,  because  the  (y,q)  =  (1,1)  term  is  much  smaller  than  the 
(y,q)  =  (1,-1)  term  in  Table  3.  The  strong  variation  in  e  during  1976-7  would 
have  continued  into  1978  if  the  great  increase  in  drag  had  not  carried  the  orbit 
away  from  resonance. 

A  second  feature  of  interest  in  Fig  4  is  the  form  of  the  variation  near 

the  central  (<J>  =  0)  resonance.  Equation  (9)  shows  that  the  variation  in 

eccentricity  is  governed  by  terms  of  the  form  sln(to  ±  <J>)  and  sln(to  +  20):  so 

cos  cos 

near  0=0  the  oscillations  should  have  the  same  period  as  to  ,  and  this  expecta¬ 
tion  is  borne  out,  as  shown  by  the  to  =  270°  values  marked,  the  first  five  of 
which  coincide  with  maxima  in  the  observational  values. 

The  most  important  feature  of  Fig  4,  however,  is  the  deficiency  in  fitting. 
This  was  a  severe  setback,  which  provoked  a  thorough  examination  of  the  modus 
operandz  of  THROE  and  PROD,  including  tests  of  omitting  alternate  points, 
running  time  backwards,  etc.  No  flaw  was  found.  So  a  reconsideration  of  the  pro¬ 
cedures  was  required.  In  the  past,  the  inclination  has  been  successfully  analysed 
for  approximately  equal  times  before  and  after  the  exact  resonance  (<J>  =  0),  and 
preferably  for  a  time  long  enough  to  include  a  few  of  the  subsidiary  oscillations 
(as  in  Fig  3).  The  eccentricity  has  in  the  past  then  been  analysed  over  the  same 
time  interval  as  i  .  The  greatest  effects  of  resonance  on  eccentricity  are  likely 
to  occur  near  the  times  when  <t>  =  <I>  and  4>  =  -  to  ,  since  the  ($  -  oj)  and  (0  +  to) 
terms  in  equation  (9)  will  then  be  nearly  constant.  Although  the  $  =  0  resonance  has 
no  relevance  for  eccentricity,  the  'eccentricity  resonances',  4>  =  +  to  ,  usually 
occur  within  a  few  months  of  the  'inclination  resonance'  and  the  analysis  can  be 
carried  through  from  before  $  =  to  to  after  $  =  -  to  (or  from  $  =  -  to  to 
>  -  '•  if  i  *  63°).  For  I97I-54A,  however,  $  —  to  occurs  in  December  1972, 

nearly  2  years  before  4  =  0;  and  4>  —  -  to  occurs  in  July  1977,  2\  years  after 

★ 

1=0  .  Thus  I971-54A  really  experienced  two  separate  'eccentricity  resonances' 

*  These  are  the  months  when  $  is  equal  to  the  mean  value  of  to  (or  -to)  over 
a  cycle:  because  of  the  variations  in  to  for  such  a  nearly  circular  orbit, 
an  exact  date  cannot  be  specified,  but  the  envelope  of  the  minor  oscillations 
in  J>  1  a)  has  its  minimum  at  about  these  times. 
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centred  more  than  4£  years  apart;  and  the  variation  in  e  near  $  =  0,  which, 
as  we  have  seen,  is  effectively  dependent  only  on  oi  ,  may  well  contain  no  use¬ 
ful  information  about  the  !5th-  and  30th-order  coefficients  in  the  geopotential. 
So  the  obvious  procedure  is  to  abandon  previous  practices  and  to  treat  the  two 
e  resonances  separately.  Values  of  e  near  4>  =  0  can  be  excluded  unless  they 
are  helpful. 

In  analysing  the  two  resonances,  not  all  the  terms  in  equation  (9)  are 
needed:  the  small  terms  will  have  a  negligible  effect  at  the  'other'  resonance, 
where  they  will  only  generate  small  periodic  terms.  Thus  for  the  first  reson¬ 
ance  we  shall  expect  to  use  the  (y,q)  =  (!,!),  (1,-1)  and  (2,1)  terms,  the  (1,-1) 
terms  being  included  because  they  are  so  large;  and  for  the  second  resonance  just 
the  (1,-1)  and  (2,-1)  terms  would  be  used.  Again,  however,  it  may  be  necessary 
to  add  (y,q)  =  (0,1)  terms. 


In  the  first  e-resonance  95  orbits  were  used  initially  (MJD  41633-42291). 
The  fitting  with  (y,q)  =  (1,1), (1,-1)  and  (2,1)  was  fairly  satisfactory,  with 
e  =  1.23;  but  there  was  a  considerable  improvement  when  the  (0,1)  terms  were 
added,  z  being  reduced  to  0.74  and  the  standard  deviation  of  the  (1,1)  terms 
being  reduced  by  40%.  The  addition  of  other  terms,  such  as  (3,1)  and  (2,-1), 
was  tested,  but  gave  no  advantage.  The  fitting  with  (y,q)  =  (1 ,!),(! »  — 1),(2, 1) 
and  (0,1)  was  quite  satisfactory,  except  that  the  last  16  values  did  not  fit 
well.  When  these  values  were  omitted,  z  decreased  from  0.74  to  0.58,  the 
values  of  the  (1,1)  coefficients  changed  by  less  than  0.4  sd,  and  the  standard 
deviations  decreased  by  more  than  10%.  This  shortened  fitting,  of  79  values, 
seems  slightly  preferable.  A  further  shortening,  to  7!  values,  was  also  tried. 
There  was  a  decrease  in  e  ,  to  0.53,  and  a  slight  decrease  in  most  but  not  all 
of  the  standard  deviations.  The  relative  merits  of  these  last  two  fittings  are 
debatable,  but  the  argument  against  the  shorter  run  (namely  the  increased  danger 
of  overfitting)  led  to  a  preference  for  the  79-value  fitting.  This  is  taken  as 
the  final  choice  and  shown  in  Fig  5. 

In  the  second  e-resonance,  117  orbits  were  included  (MJD  42697-43522).  The 
fit  with  (y,q)  =  (1,-1)  and  (2,-1)  only  was  fairly  good,  with  e  =  1.66,  but  the 
inclusion  of  the  (0,1)  terms  again  improved  the  fit,  z  being  reduced  to  0.94, 
with  large  reductions  in  two  of  the  four  standard  deviations.  This  fitting  is 
shown  in  Fig  6  and  was  taken  as  final,  since  the  addition  of  other  terras, 
such  as  (3,-1)  and  (1,1)  gave  no  advantage;  reducing  or  increasing  the  number  of 
orbits  to  97  or  131  was  also  tried,  but  found  to  be  of  no  advantage. 
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Though  the  fittings  shown  in  Figs  5  and  6  are  entirely  satisfactory,  we 
have  to  ask  (a)  why  the  (y,q)  =  (0,1)  terms  are  needed,  and  (b)  whether  their 
presence  disturbs  the  values  of  the  lumped  harmonics  obtained.  The  (0,1)  terms 
represent  an  oscillation  with  the  same  period  as  u  ,  so  the  need  for  them 
suggests  that  the  w-dependent  oscillation  in  eccentricity,  which  is  partially 
removed  from  the  US  Navy  values  of  e  ,  has  not  been  correctly  restored.  In 
some  previous  analyses  such  an  error  was  certainly  present  and  (0,1)  terms  were 
included  to  correct  it.  But  ELTRAN2  was  intended  to  solve  the  problem.  It  has 
not  done  so,  and  the  question  will  be  pursued  further.  In  answer  to  the  second 
question,  it  seems  that,  since  no  pure  w  terms  due  to  resonance  appear  in  the 
equation  (5),  there  is  no  reason  why  the  fitting  of  an  u  term  should  degrade 
the  accuracy  of  the  lumped  harmonic  coefficients  obtained:  on  the  contrary, 
the  inclusion  of  a  term  which  greatly  improves  the  fitting  should  greatly  improve 
the  values. 

If  the  ai  term  arises  from  the  failure  to  remove  oi-dependent  terms  from 
the  original  data,  the  values  found  for  the  coefficients  of  the  (0,1)  terms  from 
analysis  of  the  first  and  second  e-resonances  and  all  269  values  should  show 
similarities.  In  fact  the  values  of  the  sin  u>  and  cos  w  coefficients  from 
the  first  resonance  (Fig  5)  are 

-8.0  ±0.9  and  -10.8  +1.1; 

the  corresponding  values  from  the  second  resonance  (Fig  6)  are 

-10.2  ±  0.8  and  -8.5  ±  1.4  ; 

while  the  values  from  the  269-orbit  fitting  (Fig  4)  are 

-13.0+1.0  and  -11.8  ±  2.1  . 

.  -9 

These  three  pairs  of  values  (all  are  x  10  )  agree  with  each  other  to  within 

2.6  sd  (sin)  and  1.0  sd  (cos),  and  so  can  be  regarded  as  consistent. 

6 . 3  Values  of  the  lumped  harmonics 

The  values  of  the  lumped  harmonics  obtained  from  the  three  fittings  of  e 
are  given  in  Table  4.  The  preferred  pairs  of  values  are  underlined. 

Table  4 

Values  of  lumped  harmonics  obtained  from  fitting  of  e 


Uccrhi  (l.-l)  (2,-1) 

r*totuuic»  (0,1) 

(117  valu«») 


I 


The  results  in  the  first  row,  from  all  269  values,  have  already  been  discounted 
because  of  the  poor  fit;  but  some  of  the  lumped  harmonic  values  are  in  agreement 

-  -  -  ]  2 

with  those  from  the  two  divided  fittings.  The  preferred  values  of  (C,S)jj*  are 
those  from  the  second  resonance,  of  course,  because  this  is  the  main  term  at  that 
resonance.  The  values  of  (C,S)j^  were  expected  to  be  -  and  are  -  of  poor  accu¬ 
racy,  because  their  influence  on  the  resonance  was  only  about  1/20  of  that  of  the 
(C,S)  terms  (Table  3). 

On  taking  the  mean  of  each  underlined  pair  of  values  of  |c|  and  |s|  in 
Table  4,  and  multiplying  them  by  the  appropriate  F  factors,  etc,  in  equation  (9), 
we  find  that  the  relative  magnitudes  of  the  four  harmonics  are  9:100:13:9,  as  com¬ 
pared  with  the  a  priori  order-of-magnitude  expectation  (Table  3)  of  5:100:8:12. 

Thus  the  order  of  magnitude  estimates  were  excellent,  being  within  a  factor  of  2 
of  the  real  values. 

7  THE  LUMPED  HARMONICS  EXPRESSED  IN  TERMS  OF  THE  INDIVIDUAL  C„  ,  S„ 

_ Jim  Jim 

,q»k  q,k  .  . 

The  lumped  harmonics  Cm  and  Sjj,  are  linear  functions  of  the  indi¬ 
vidual  harmonic  coefficients  C and  S£m  given  by  equations  (6).  The  numeri- 
cal  values  of  the  coefficients  Q^’  in  equations  (6)  for  1971-54A  have  been 
evaluated  with  the  aid  of  the  computer  program  PROF,  to  express  each  of  the  6 
pairs  of  lumped  harmonics,  given  in  Tables  2  and  4,  in  terms  of  the  individual 
harmonics.  The  numerical  values  of  the  Q^’  decrease  when  the  degree  Jl  is 
increased  far  enough,  and  only  a  limited  number  of  terms  need  to  be  evaluated, 
depending  on  the  accuracy  of  the  values  of  the  lumped  harmonics.  Here  we  have 
ignored  terms  which  are  not  expected  to  contribute  (relative  to  the  main  term)  an 
amount  more  than  half  the  proportional  error  in  the  value  of  the  more  accurate  of 

the  two  lumped  coefficients.  Thus,  for  the  first  pair  of  lumped  coefficients,  the 
_0, I  _9 

more  accurate  is  C  ,  which  has  a  value  of  (-16.4  ±  0.24)  x  10  ,  so  that  the 

^  0  1- 

proportional  error  is  0.015.  The  expected  value  of  Q  C  relative  to  the 

0  1  2  x,  Jim  _ 

first  term  (which  ;s  the  largest)  is  Q  ’  (JIq/JI)  ,  if  we  assume  is  of 

—52  .  01 

order  10  /?“  .  Thus,  since  l  =  15  for  ,  the  series  is  truncated  when 

0  i  t  0  15 

Q„  (l5/»’.)“  ■  0.0075  :  this  applies  for  £  *  33  ,  so  the  series  is  taken  to  Jl  =  31. 

Under  these  criteria  the  equations  for  the  12  lumped  harmonics  are  as  follows. 

C  :  r",5,l5  +  °-5Mr:,7,15  +  °-  327S9,l5  +  °-22l521,15  +  °-,54623,l5 

+  °-l,0S5,15  +  °-079r:27,i5  +  °-058'C29,.5  +  °-°42Sl,15 

=  (-16.4  ±  0.24)  x  10-9.  (10) 

_0,l 

:  the  equation  is  the  same  as  (10),  with  S  instead  of  C  ,  and  the 

-9 

numerical  value  on  the  right-hand  side  changed  to  (-5.4  ±  0.15)  x  10  .  (11) 


v  . 


C30,30  +  °-430C32,30  +  °-213C34,30  +  °-100C36,30  =  ("8-2  * °-6>* 10  •  ('2> 

the  equation  is  the  same  as  (12),  with  S  instead  of  C  ,  and  the 

•  _Q 

numerical  value  on  the  right-hand  side  changed  to  (11.1  ±0.8)  x  10 


'16,15 


+  ,*I77C18,15  + 

'•250C20,15  +  *■ 

•270C22,15  + 

'•256C24,,5 

+  '•220S6,15  + 

'•,69528.I5  +  ' 

•108So,15  + 

,.040C32j,5 

+  °-970S4,«5  + 

II 

IT) 

vO 

rn 

ICJ 

o> 

O' 

oo 

o 

(-  92  ±  48) 

i 

o 

X 

the  equation  is  the  same  as  (14),  with  S  instead  of  C  ,  and  the 
numerical  value  on  the  right-hand  side  changed  to 
(-  170  ±  56)  x  io-9 


C 1 6 , 15  +  °-957C18,15  +  °-839C20,15  +  °-714C22,.5  +  °-598C24,15 

+  °-4%526,.5  +  °-4!0S8,15  +  °-337So,15  +  0-276S2,15 

+  °-226S4,15  +  °-185S6,15  +  °-,51S8,15  +  °'123S40,15 

=  (-  62.9  ±  2.6)  x  IO"9 

the  equation  is  the  same  as  (16),  with  S  instead  of  C  ,  and  the 
numerical  value  on  the  right-hand  side  changed  to 
(-  53.4  ±  1.6)  *  IO-9 

Si, 30  *  l-040S3,30  *  °'988S5,30  *  °-90,S7,30  *  0'8I7S9,30 

*  0.728C41;30  *  0.644C43>30  4  0.5670^  4  O.496C47>30 

4  0.433C49  m  -  <15  !  8)  *  IO'9 

the  equation  is  the  same  as  (18),  with  S  instead  of  C  ,  and  the 
numerical  value  on  the  right-hand  side  changed  to 
(-  39  ±  12)  x  10~9 


C  3 1 , 30  +  °*874C33,30  +  °'695C35,30  +  °*531C37,30  +  °‘394C39,30 

=  (-  8  ±  3)  x  10~9 


:  the  equation  is  the  same  as  (20),  with  S  instead  of  C  ,  and  the 
numerical  value  on  the  right-hand  side  changed  to 

(-  13  ±  7)  x  io'9.  (21) 

Equations  (10)  to  (21)  provide  12  equations  of  condition  for  the  individual 
coefficients  of  order  15  and  30,  which  will  be  used  in  future  determinations  of 
the  individual  coefficients. 


qq  n 


19 


8  TESTS  OF  COMPREHENSIVE  GRAVITY  FIELD  MODELS  BY  COMPARISON  OF  THEIR  LUMPED 
HARMONICS  WITH  OBSERVATIONAL  VALUES  FROM  1971-54A 

In  recent  years  a  number  of  gravity  field  models,  complete  up  to  a  certain 
order  and  degree,  have  been  determined  from  a  combination  of  satellite  observa¬ 
tions,  terrestrial  gravity  measurements  and  altimeter  data.  Of  these  models  the 

26 

Smithsonian  Standard  Earth  IV. 3  goes  to  order  and  degree  24,  the  European 

27  1 

model  GRIM  2  goes  to  order  and  degree  30,  and  the  Goddard  Earth  Model  10B 

2 

goes  to  order  and  degree  36;  while  GEM  1 0C  consists  of  GEM  1  OB  up  to  order  and 
degree  36,  together  with  31000  extra  coefficients  determined  from  the  altimeter 
geoid  and  going  up  to  order  and  degree  180. 

Inspection  of  equations  (10)  to  (21)  shows  that  terms  of  degree  higher 

than  24  are  involved  in  all  the  equations,  so  comparisons  with  SSE  IV. 3  cannot 

be  pursued.  Comparisons  with  GRIM  2  are  also  fruitless,  because  the  15th-order 

3  4 

terms  in  GRIM  2  rely  heavily  on  our  1975  sets  of  individual  harmonics  ’  which 
utilize  our  analysis  of  197I-54A  for  the  years  1972-4.  In  fact  the  values  of 

q_  0 ,  1  n  0  ,  1 

10  C  and  10  from  GRIM  2  are  -19  and  -4,  as  compared  with  -16.4  +  0.2 

and  -5.4  +0.1  in  Table  3.  The  agreement  is  within  the  expected  standard 
deviation  of  GRIM  2. 

This  leaves  only  GEM  10B  (or  GEM  10C  if  terms  beyond  P.  =  36  are  required) 

for  comparison.  In  making  the  comparisons,  the  last  four  equations  (18)— (21), 

will  be  omitted,  because  the  lumped  values  are  not  accurate  enough  to  provide  a 

satisfactory  test  -  their  standard  deviations  range  between  30%  and  55%.  With 

these  omissions,  only  one  of  the  pairs  of  equations,  namely  (16)  and  (17), 

includes  terms  of  degree  higher  than  36,  and  these  are  small,  so  that  the 

180  x  180  field  GEM  10C  is  not  really  being  tested  and  the  comparison  is  with 

GEM  10B.  Substituting  the  values  of  C.  ,  S„  from  GEM  1 0B  into  equations  (10), 

Urn  Jim 

(II),  (14),  (15),  (16)  and  (17),  ignoring  the  small  terms  involving  GEM  10C  in 
(16)  and  (17),  we  obtain  values  for  the  15th-order  lumped  coefficients  which  are 
compared  with  those  from  the  analysis  of  197I-54A  in  Table  5.  It  was  tentatively 

suggested  in  Ref  28  that  the  individual  coefficients  in  GEM  10B  might  be  given 

-9 

standard  deviations  of  5  x  10  ,  and  the  standard  deviation  of  the  lumped 

coefficients  has  been  estimated  on  this  basis. 

Table  5  shows  that  the  second  pair  of  coefficients  (C,S)j’^  is,  nominally, 
better  determined  from  GEM  I0B  than  from  the  resonance  analysis.  This  was  to  be 
expected,  because  the  (y,q)  =  (1,1)  terms  have  very  little  influence  on  the  reson¬ 
ance,  as  indicated  by  Table  3.  So  the  second  pair  of  lumped  values  from  I97I-54A 
is  of  no  use  as  a  test  of  GEM  I 0B. 


i 
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Table  5 


Values  of  15th-order  lumped  harmonics  from  1971-54A  resonance, 
compared  with  corresponding  values  from  GEM  10B 


q  0,1 

109C,5 

o  0.1 

10  5, 5 

Q  1  ,0 

l0C,5 

q  1.0 

.0  s15 

9-_1  ,2 
10  C,  5 

9-_1  ’2 
10  S 1  5 

1971-54A 

-16.4  ±  0.2 

-5.4  ±  0. 1 

-92  ±  48 

-170  ±  56 

-63  ±  3 

-53  ±  2 

GEM  10B 

-22  ±  6 

-8  ±  6 

-77  ±  18 

-55  ±  18 

-64  ±  10 

-46  ±  10 
_ ! 

For  the  other  four  coefficients  in  Table  5,  Llie  resonance  values  are  nomin¬ 
ally  much  more  accurate  than  GEM  10B,  and  provide  the  required  test.  GEM  10B 
emerges  from  the  comparison  extremely  well,  all  the  GEM  10B  values  being  within 
1  standard  deviation  of  the  values  from  resonance  -  the  actual  differences  being 

0.9,  0.4,  0.1  and  0.7  of  the  estimated  GEM  10B  standard  deviations.  These  four 

-9 

values  have  an  rms  of  0.6,  thus  indicating  that  a  value  of  3  *  10  ,  rather  than 
-9 

5  x  10  as  assumed,  is  appropriate  as  an  approximate  standard  deviation  for 
the  15th-order  coefficients  in  GEM  10B. 

Only  one  of  the  pairs  of  30th-order  coefficients  from  the  resonance  analy- 

-  -  o  2 

sis,  namely  »  is  determined  accurately  enough  to  provide  a  real  test  of 

the  30th  order  coefficients  in  GEM  10B,  and  fortunately  the  numerical  coefficients 
in  equation  (12)  decrease  rapidly,  so  that  terms  of  degree  higher  than  36  arc  not 

needed.  The  comparisons  are  made  in  Table  6,  where  the  assumed  standard  deviation 

....  .  .  -9 

of  the  individual  GEM  I0B  coefficients  has  been  reduced  to  3  *  10  ,  as  suggested 

in  the  previous  paragraph. 

Table  6 

-  -  0  2 

Values  of  (C,S)^q  from  1971-34A  resonance,  compared  with 
corresponding  values  from  GEM  10B 


q  0,2 

q  0,2 

10 

109S„„ 

30 

30 

I97I-54A 

-8.2  ±  0.6 

11.1  ±0.8 

GEM  10B 

-8.4  ±  3.3 

11.2  +  3.3 

The  agreement  between  these  quite  independent  results  from  the  resonance 
analysis  and  GEM  10B  is  amazing,  and  should  be  treated  as  yet  another  example  of 
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very  low  probabilities  occurring  more  often  than  they  ought  to*.  If  we  "season 

our  admiration  for  a  while",  we  may  conclude  from  Table  6  that  GEM  10B  has 

passed  a  stringent  examination  with  distinction,  and  that  the  standard  deviation 

-9 

of  the  individual  values  may  be  less  than  3  *  10 

9  FURTHER  WORK 

Large  numbers  of  observations  of  I971-54A  have  been  made  by  a  variety  of 

sensors,  including  the  Baker-Nunn  cameras,  the  Hewitt  camera  at  Malvern,  the 

US  Navy  Navspasur  system  and  visual  observers.  If  orbits  were  determined  from 

all  these  observations  for  dates  near  the  two  e-resonances,  a  better  analysis  of 

the  variations  in  eccentricity  should  be  possible,  without  the  need  for  the 

(y,q)  =  (0,1)  terms.  It  would  then  be  possible  to  make  a  combined  fitting  of 
.  .  23 

i  and  e  together  using  the  SIMRES  program  ;  this  is  not  possible  at  present 
because  the  presence  of  the  (0,1)  terms  would  spoil  the  fitting  of  i  . 

10  CONCLUSIONS 

The  variations  in  the  inclination  and  eccentricity  of  1971-54A  between 
November  1972  and  January  1978,  as  given  by  the  269  US  Navy  weekly  orbits,  have 
been  analysed  to  determine  12  lumped  geopotential  coefficients  of  order  15  and 
30,  given  in  Tables  2  and  4. 

The  fitting  of  the  variations  in  inclination  is  nearly  perfect  (Fig  3)  and 

-0, 1 

yields  the  most  accurate  values  yet  obtained  for  lumped  coefficients,  C,c 

_0,2  15 

having  an  accuracy  of  1.5%,  and  an  accuracy  of  7%. 

Difficulties  were  encountered  in  fitting  the  variations  in  eccentricity 
over  the  full  5{  years,  and  the  problem  was  solved  by  separate  analyses  of  the 
two  eccentricity  resonances,  which  are  4  years  apart  (Figs  5  and  6). 

Each  of  the  12  lumped  coefficients  obtained  is  a  linear  function  of  the 
individual  harmonic  coefficients,  and  the  12  equations,  given  as  equations 
(10)  —  (21  ),  will  shortly  be  used  in  a  new  solution  for  the  individual  harmonic 
coefficients  of  15th  and  30th  order. 

Comparison  of  the  lumped  harmonics  of  15th  order  with  those  obtained  from 

the  Goddard  Earth  Model  I 0B  shows  excellent  agreement,  and  suggests  that  the 

accuracy  of  the  high-order  coefficients  in  GEM  10B  is  better  than  had  been 

-9 

expected,  being  about  3  *  10 

*  The  values  we  obtained  in  1975  for  the  lumped  coefficients  from  197I-54A  were 
(-10.3  t  1.5)  and  (15.3  ±  1.3);  so  the  cynical  suggestion  that  our  1975  values 
somehow  found  their  way  into  GEM  10B  can  be  withdrawn  before  it  is  voiced. 
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Only  one  pair  of  30th-order  coefficients  is  accurate  enough  to  provide  a 
test  of  GEM  10B.  The  values  agree  astonishingly  well,  thus  suggesting  fv’e 
possibility  that  the  GEM  IOB  coefficients  may  be  even  more  accurate  than  indi¬ 
cated  in  the  previous  paragraph. 
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variation  of  inclination  near  15th -order 


Fig  4  1971-54A:  variation  of  eccentricity  from  IMov  1972  to  Jan  1978 
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Fig  5  1971-54A:  variation  of  eccentricity  near  the  (y,q)  =  (1,1)  resonance 
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variation  of  eccentricity  near  the  (7,q)  -  (1,— 1)  resonance 


